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Abstract 

] The photon emission rates from the quark gluon plasma have been studied considering LPM 

suppression effects. The integral equation for the transverse vector function (f (p_i_)) that consists 
of multiple scattering effects has been solved using self-consistent iterations method. Empirical 
fits to the peak positions of the distributions from iteration method have been obtained for 
bremsstrahlung and aws processes. The variational approach for f(p_i_) calculation has been 
. simplified considerably making some assumptions. Using this method, the photon emission rates 

at finite baryon density have been estimated. The LPM suppression factors for bremsstrahlung 
and aws processes have been obtained as a function of photon energy and baryon density. The 
effect of baryon density has been shown to be rather weak and the suppression factors are similar 
to the zero density case. The suppression factors for aws processes can be taken at zero density, 
| whereas the bremsstrahlung suppression can be taken at zero density multiplied by a density 

O ■ dependent factor. 
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Electromagnetic processes such as photons and dileptons production are considered to be impor- 
q^' tant signals of formation of the quark gluon plasma (QGP) in the relativistic heavy ion collisions. 

The experimental data by WA98 collaboration ^ of direct photons from Pb+Pb collisions at 
CERN SPS is of paramount importance which boosted interest in the static emission rates and 
the convoluted photon yield calculations. A comprehensive review on the theoretical and exper- 
imental studies of the photon emission has been given in Q. The photon emission rates from 
hot hadron gas, the QGP phase and also the prompt photon emission have been investigated 
in great detail (see references in Q). The processes contributing to photon emission from QGP 
phase at (HTL) effective one loop level are the quark-antiquark annihilation into a photon and 
a gluon and the absorption of a gluon by a quark (anti quark) emitting photon (see Fig.l of 
Q ) . The processes of bremsstrahlung and aws that arise at effective two loop level contribute 
at the same order as one loop processes Q. These processes (see Fig. 5 in Q) contribute at 
the leading order 0(aa s ) owing to the collinear singularity that is regularized by the effective 
thermal masses. The higher order multiple scatterings however can not be ignored as these may 
also contribute at the same order as the one and two loop processes Q || ](| [j], ^ [|. Further, 
multiple soft scatterings of the fermion during photon emission reduce the photon coherence 
lengths, known as Landau-Pomeranchuk-Migdal (LPM) effect. The photon emission rates are 
suppressed owing to the LPM effects [§, @, §]• Using the hydrodynamical models for the 
expansion of the plasma and convoluting the expansion with the photon emission rates, one 
obtains the photon transverse momentum spectrum. These space time evolution of the plasma 
through QGP, mixed and hadron phases is complicated if the plasma is formed at finite baryon 
density. These physical conditions of the plasma affect both the basic photon emission rate and 
the evolution of the plasma and thus the space time integrated photon spectrum. In view of this, 
the formalism of Aurenche et. al, of photon emission rates at two loop level has been recently 
extended to a chemically unsaturated plasma at finite baryon density [Gi]] . These emission rates 
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(meaning same as production rates in this work) and the space time (3+1 dimensional) expan- 
sion dynamics of the plasma at finite baryon density have been recently reported |Tl| , |l2"fl . The 
photon spectra have been shown for various cases of plasma at finite baryon density for SPS 
and RHIC energies and the baryon free case in @ included the LPM effects. 



The photon production rates from bremsstrahlung and the aws processes have been estimated 
by Aurenche et. al, [||. In a remarkable simplification, these rates are expressed in terms of 
simple one dimensional momentum integrals and the dimensionless quantities Jt,Jl H- These 
values are not very sensitive to baryo chemical potential and they weakly depend on the thermal 
masses for gluon and fermions through m 2 /m 2 a |l(J. The differential photon emission rate per 
unit volume of energy ko is given by (denoted by 1Z° for simplicity of later use) , 

A« = ^^n B (k )^(J T - Jl) hM (1) 
I b (k)=J^ dp(p 2 + (p + k) 2 ^[(n f (p) -n f (p + k)) + (n f -n f (p + k)} (2) 
I a (k) = f°° dp (p 2 + (k - p) 2 ) [1 - n f (p) -n f (k- p)\ (3) 



In above, the subscripts (b, a) denote the bremsstrahlung and the aws processes and we consider 
a two flavor three color case with a s = 0.2. These rates in Eq. 1 do not include LPM suppression 
effect. The effects of finite baryon density on the photon emission rates have been discussed 



in [11] and these results are summerised as follows. The bremsstrahlung radiation is affected 
whereas the aws process is insensitive to the baryon density presence. The bremsstrahlung 
radiation from quark is enhanced and from anti quark is suppressed at finite baryon density. As 
mentioned earlier, it has been shown that multiple scatterings during the photon emission for 
the bremsstrahlung and aws processes cannot be ignored. The exchange of a soft gluon of finite 
thermal mass and the collinear singularity (arising when photon is emitted parallel to quark 
momentum) regularized by the thermal masses are essentially the factors that make the two 
loop processes contribute at the same order as the one loop processes. It has been realized []|, [?j 
that for similar reasons, multiple gluon exchange processes in terms of certain ladder diagrams 
also contribute to the same order. The diagrammatic analysis of all the processes contributing 
to the same order has been performed and these contributions have been summed in (7j. This 
leads to suppression of photon emission rates as shown by the LPM suppression factors in Fig. 
7. of . One can also obtain the suppression of the total emission rate by using the empirical 
expressions in Eq.(l.lO) of H and the two loop rates of Eq. 1 above without LPM effects. It has 
been shown that the photon radiation from bremsstrahlung is strongly suppressed to ~ 0.2 at 
very low k/T values. The suppression becomes weaker at higher photon energies and approaches 
~ 0.86 for k/T > 10. In contrast, the aws is unaffected at low k/T values, but falls almost 
linearly up to k/T ~ 10. The bremsstrahlung from anti quark is same as that of quark for a 
baryon free case. In the presence of finite baryon density, the bremsstrahlung from quark will be 
different from the anti quark. Therefore, in order to compare with experimental photon spectra, 
it is necessary to understand how the suppression factors change with baryon density. The 
chemical potential dependence of the emission rates are contained in the population functions 
of the quarks (anti quarks). Further, for comparing with experimental data it is necessary to 
determine these suppression factors (beyond k/T ~ 10) up to k/T ~ 20 i.e., k ~ 5 GeV for 
T=0.25GeV case. The differential photon emission rate (denoted by 1Z) that includes the LPM 
effects is given by (for details see ||), 
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n f {k+p\\){l-n f {p\\)) 2p ± -Kf(pj 



(4) 



In the above k = m^/mj^ (=2 for fi = case) and the subscripts (b, a) are for 
bremsstrahlung and aws as in Eq. 1, with different kinematic domains and appropriate dis- 
tribution functions. The value of k is baryon density dependent and is determined by evaluating 
these thermal masses as listed in the Table as a function of /i q . In the above equation, 5ftf(pj_) 
is the real part of a transverse vector (amplitude) function which consists of the LPM effects 
due to multiple scatterings. This can be taken as transverse vector (p_i_) times a scalar function 
of transverse momentum p± . The sign ~ denotes the dimensionless quantities in units of Debye 
mass mo as defined in ||]. The function pj_ • 5ftf (pj_) is determined by the collision kernels 
(C(q_|_)) in terms of the following integral equation. 



2p ± = iSE(p±,p\\,k)f(p±,p\\,k) + 



(2vr)2 



f(p±iPiijfc)-f(p± + q±,Pii,AO ^(q 



(5) 



<?(qj 



+ 



(1 - (<?/q) 



SE(p±,p h k) 



kT 



2p\\(k + p\\) 



p\ + K 



(6) 
(7) 



In above, the T1l,T1t are the longitudinal and transverse parts of (dimensionless) thermal self 
energies of the gauge fields and 5E{p±,p», k) is the energy difference between the relevant 
states of the system before and after the photon emission ||]. Based on the sum rules for the 
thermal gluon spectral functions, Aurenche, Gelis and Zaraket obtained an analytical form for 
the collision kernel as given by Eq. 44 of Q. This enormously simplifies the photon emission 
rate calculations, as it circumvents the need for evaluating the integral in Eq. 6 to obtain the 
collision kernel. In the present work, we adopt this simple analytical form for the collision kernel, 
with appropriate factors properly taken care of. We have computed the integral in Eq. 6 for 
the collision kernel and compared with the analytical form given in ||. We found these to be 
agreeing, except near = where the collision kernel of Eq. 6 falls to zero and can be fitted to 
another analytical form for use in calculations. This result has been surprising as the analytical 
form in || is based on very general assumptions satisfied for gluon thermal (HTL) self energies 
and sum rules for its spectral functions. The analytical form of || for the full integration region 
can also be used in the calculations. The divergence of kernel in || at q_i_=0 is compensated 
by the vector function difference and the d 2 q_|_ together. Further, use of rationalized momenta 
in terms of debye mass cancel the 3m 2 , terms in the analytical form given by Aurenche et. al.,. 
Therefore, the C(qj_) function just scales as k times the analytical forms which are independent 
of chemical potentials. 

The function pj_ • 3?f(pj_) as a function of |pj_| has to be solved self consistently for each set of 
{pii , k, k, T} values. In the present work, we have solved the equation by iterations at a fixed 
temperature of T=0.25GeV and obtained 2p^ • 5Rf(pj_) distribution. For small pn and k values 
the iterations converge very fast. For small k and large P\\(> 2.5) the convergence is slow. The 
peak positions of these distributions have been obtained for each set of {pu,k,tz}. It has been 
noticed that the peak positions of iteration method (Aj) can be very well approximated by, 
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The value of (3 parameter is, j3 = —0.32 for bremsstrahlung and aws away from the py limits. 
(The value of (3 = —0.16 for aws near the py limits has been a good approximation.) Further, a 
lowest cut-off value of A™ 171 = 0.32 was taken, which is necessary near py limits or very small k 

values, especially for aws case. The f distributions have been studied for various values of quark 
chemical potentials ranging from 0-2GeV. It has been found that the peak positions are not very 
sensitive to the chemical potentials. The dependence of the chemical potential is contained in 
k factor of the integral equation which weakly depends on as shown on the table. Figure 
l(a, b, c, d) show the distributions of p± • (p_i_) for typical cases of {fc,py , k} values. The Figs. 
l(a,b) are for bremsstrahlung and Figs. l(c,d) are for aws with the parameter values shown in 
figures. The black curves show the results from iterations for zero baryo chemical potential and 
the blue curves show for the case of /i = l.OGeV. As shown in these figures, the A{ values do 
not change significantly with p. The values are generally less than unity for most parameter 
values for aws case. The empiricism has been derived from a study of several plots of the 
amplitude function for various parameters and at various stages during iterations. The resulting 
2pj_ • 5if(p_i_) distribution from the iteration method can be used to perform the photon rate 
calculation using Eq.(4), though the method is not efficient. Results of this method for low 
photon momentum and high py values become unreliable when the convergence is poor, though 
the convergence can be improved by standard methods. However, higher py values (py > 10T) 
contribute lesser to the emission rate owing to the population functions. Therefore, in the 
following we adopt a different method for the photon emission rates. 

We have followed the variational approach discussed in detail in Eqs. [4.20-4.29] of || for solv- 
ing the integral equation. We found that this method is remarkably simple to solve for the 
2pj_ ■ 5ftf (pj_) distributions, though highly intensive in computing time. The choice of trial func- 
tions of H] do not allow the integrals to be analytically tractable and hence can be computed 
only numerically. Further, proper choice of the scale constant of the variational method A v is 
necessary for use in these trial functions. The A v value should be around the peak position of 
the 2pj_ • Kf(pj_) distributions. For this purpose, a constraint integral has been given in Eq. 
4.24 of H . In the present work, we simply this approach as follows. 

(i) The scale constants (A v ) are taken from empirical expression of Eq. 8, i.e., A v (p\\,k) = A\ 
for brensstrahlung , A v (p»,k) = Ai for aws. 

(ii) These A4 values are already shown to be weakly dependent on baryon density. Therefore 
the scale constants are taken same as for the \x = case. 

(iii) The collision kernel of Eq. (6) has been replaced by the analytical forms as mentioned 
in the last section. 

These three considerations grossly simplify the problem and bring out the full advantage of 
the variational method suggested in ||. We have calculated the 2pj_ • 5?/(pj_) distributions for 
various sets of {pn , k, k}. We compared these with the iteration method mentioned earlier. These 
distributions are shown by red colored curves in Figs, l(a-d) using the dimension (N r ) 12, of 
the set of trial functions. The pink curve is the corresponding distributions for \i = 1 GeV case. 
As shown in the figures, the agreement between these two methods is good. This agreement 
remains valid even for low value of N r . At finite baryon density, use of dimension N r = 5 or 6 
is just sufficient for the photon emission rate calculations using A v (pu,k). Further, for a choice 
of N r =12 or 15, these distributions will not be sensitive to the A v values used. This is expected 
in any basis expansion method for sufficiently large dimension of the model space. Therefore, 



prudent choice of set of functions and parameters seems necessary to get converging results only 
for lesser dimensions that render computations easy. This has been verified by comparing the 
distributions for N r =4 using Ao(pii, k) for a few {pn, k}, with the results for N r =12 or 15 using 
any fixed value of A v . It has been observed that use of any constant value of A v with N r =8 
also gives acceptable results, provided the constant A v value is neither too large nor small as 
compared to the peak positions. It can be larger by a factor two and can be smaller by a factor 
of half than the value given by the prescription for k). A priori, the A v value for use in 

calculations for any arbitrary case is not known and therefore the aforesaid prescription serves 
as a general guideline. 



The photon emission rates have been calculated for bremsstrahlung and aws using the variational 
approach for various values of /x , as denoted by TZ^, 7Z a in Eq. 4. The emission rates for effective 
two loop processes without LPM effects have also been calculated using Eq. 1 and are denoted 
by nl,K° a . We obtained the suppression factors as defined by the ratio of the rates from these 
two methods. However for each \x value, we normalised the suppression factors to unity for the 
aws process at very low k value ( k ~ 0.2T). The bremsstrahlung suppression versus k for all jj, 
values are asymptotically normalised to zero density case, i.e., at k ~ 20T. These normalization 
factors are obtained from /&(//), f a {^) defined in the following equations. The normalised (to 
zero density) suppression factors (S&, S a ) as functions of k,jx are thus defined as, 



Sb(k^) = [^) with m = \^\ (9) 

=20T 





=0.2T 

These normalised suppression factors are shown in Figs. 2(a,b), for each \i (quark chemical 
potential) as captioned in the figures. The effects of baryon density on these suppression factors 
is weak only for low baryon density and the qualitatively these are very similar to the results 
shown in |§. Our suppression factors agree quite well with the results of (§| for baryon free case. 
The normalization factors /&, f a are shown in the table as a function of quark chemical potential. 
The values of various thermal masses and k values are also shown in the table. As seen in the 
table, the normalization factors f a for aws are very close to 1.0 showing that the baryon density 
has no effect. The aws suppression factors in Fig. 2(b) also show only weak dependence on 
baryon density. Therefore in total the aws process is insensitive to baryon density and one 
can use the suppression factors for zero density case for use in photon spectrum calculations. 
However, for the bremsstrahlung case the density effect is not totally weak. It can be seen from 
the table that the /&(//) value varies considerably for large chemical potential values. The total 
suppression curves as a function of k for the bremsstrahlung and aws at a given jj, are thus given 
by multiplying with ^'^^ the respective curves in Fig. 2(a,b) to undo the normalization. 
For low density, one can use the suppression factors as a function of k for zero density case 
multiplied by a density dependent factors f^ b '^=o) • ^ may be of interest to know from the 
table that nf^fi) is roughly a constant (=0.214) and therefore the density dependent factor 
can be chosen as kq/k. These calculations can be extended to an important case of the plasma 
conditions at RHIC collision, where the k values can be very different from the values in present 
study. 
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Conclusion 



The photon emission rates from the quark gluon plasma have been studied considering LPM 
effects. Self-consistent iterations method has been used to solve the integral equation for the 
3?f (pj_) distributions. The peak positions of these distributions have been fitted by an empirical 
expression for various parameter sets. These peak positions are observed to be rather insensitive 
to the baryon density. The variational approach has been adopted with some simplifications. 
Analytical forms for collision kernel have been used in both the iterations and variational meth- 
ods. The scale constants have been used from empirical expressions from iteration method. The 
suppression factors for bremsstrahlung and aws have been obtained as a function of photon en- 
ergy at different baryon densities. The normalised suppression factors are qualitatively similar 
to the zero baryon density case. The normalization factors show the density dependence for 
bremsstrahlung. The aws process did not show significant density dependence. 
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Figure 1: The 2pj_ • 3ftf(p_|_) distributions from iteration method (IM) for bremsstrahlung 
(Figs.(a,b)) and aws (Figs.(c,d)) processes. The black curves are for the case of zero density and 
the blue curves are for finite density with 1.0 GeV quark chemical potential. The corresponding 
cases from variational method (VM) are shown in red and pink color. The pit, k values (in GeV) 
used are shown in figures. 8 




Figure 2: (a) bremsstrahlung and (b) aws suppression factors versus photon energy. The sup- 
pression factors have been normalised to zero density case asymptotically for bremsstrahlung 
and for aws the normalization is at k = 0.2T. The true curves at nonzero density can be 
obtained from the normalization factors given in table, as discussed in text. 
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